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Exact particle- like static, spherically and/or cylindrically symmetric solutions to the equations of 
interacting scalar and electromagnetic field system have been obtained. We considered Freedman- 
Robertson- Walker (FRW) space-time as an external homogenous and isotropic gravitational field 
whereas the homogeneous and anisotropic Universe is given by the Godel model. Beside the usual 
solitonian solutions some special regular solutions know as droplets, anti-droplets and bags (confined 
in finite interval and having trivial value beyond it) have been obtained. It has been shown that in 
00 ' FRW space-time equations with different interaction terms may have stable solutions while within 

Q\ , the scope of Godel model only the droplet-like and the hat-like configurations may be stable, if they 

■ are located in the region where g 00 > 0. 
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1. INTRODUCTION 



The concept of soliton as regular localized stable solutions of nonlinear differential equations is being widely utilized 
in pure science Q. One of the fields to apply the soliton concept is the elementary particle physics, where the soliton 
solutions of nonlinear field equations are used as the simplest models of extended particles • Development of general 
| relativity (GR) and quantum field theory (QFT) leads to the increasing interest to study the role of gravitational 
y—{ • field in elementary particle physics. To obtain and study the properties of regular localized solutions to the nonlinear 
t-H I classical field equations is motivated mainly by a hope to create a consistent, divergence-free theory These solutions, 
OO , as was remarked by Rajaraman Q give us one of the ways of modeling elementary particles as extended objects 
' with complicated spatial structure. In such attempts it is natural to treat the field nonlinearity not only as a tool 
O for avoiding the theoretical difficulties (such as singularities) but also as a reflection of real properties of physical 
O" 1 system. It should be also emphasized that the complete description of elementary particles with all their physical 
, characteristics (e.g., magnetic momentum) can be given only in the framework of interacting field theory ^j. Let 
us remark that the choice the field equations is one of the principle problems in nonlinear theory. At present there 
is no criterion for the selection of interaction Lagrangian and any Lorentz-invariant combination of field functions 
can be considered as such Lagrangian. The purpose of this paper is to present some new results (anti-droplets and 
hats), in addition to those illustrated in Q for an interacting system of scalar and electromagnetic fields, confining 
ourselves to static, spherically and/or cylindrically symmetric configurations since the effect of gravitational fields on 
the properties of regular localized solutions significantly depends on the symmetry of the system. 
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2. FUNDAMENTAL EQUATIONS 

So that the field equations possess regular solutions it is necessary to introduce nonlinear terms, describing the 
field interactions, in the Lagrangian. We consider the nonlinear generalization of the theory that is related to the 
introduction of direct interaction between neutral scalar and electromagnetic fields. The decay process like 7r° — ► 27, 
described by the effective Lagrangian @ 

Lint = <Pir n Faf3F* a P , 

indicates to the possibility of such generalization. Thus we consider a system with the Lagrangian 

L = \<P, a ?> a - -^F a pF a ^(ip) (2.1) 

Z 1D7T 

where ^(y) is some arbitrary function characterizing interaction between the scalar (tp) and electromagnetic {F^ u ) 
fields takes the form 

#(<p) = 1 + \$(<p). 
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As is seen, for A = 0, ^(<p) = 1 and we have the system with minimal coupling. Note that the Lagrangian ( |2.l| ) 
describes the system of fields with positive definite energy if ^(tp) > 0. This kind of interaction has been thoroughly 
discussed in ||. Let us write the scalar and the electromagnetic field equations corresponding to the Lagrangian (pd|) 
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Q/ % = 0, 
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(2.2) 



-g dx v 



The corresponding energy-momentum tensor reads 
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<P,n<P' 



-gF^(tp) = 
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F aP F a ^{y) 



(2.3) 



(2.4) 



3. SPHERICALLY SYMMETRIC CONFIGURATIONS 



A. Solutions in FRW Universe 



As an external homogenous and isotropic gravitational field we choose the FRW space-time. This Universe is very 
important as the corresponding cosmological models coincides with observation. The interval in the FRW Universe 
in general takes the form ||,|l(J 



ds 2 



dt 2 - R 2 {t) 



dr 2 



1 — kr 2 



■{dti 



2 1 sin 2 



(3.1) 



Here R(t) defines the size of the Universe, and k takes the values and ±1. We consider the simple most case putting 
R(t) = R = constant, which corresponds to the static FRW Universe. In static case k = corresponds to usual 
Minkowski space, k = +1 describes the close Einstein Universe Jll[] and k = — 1 corresponds to the space-time with 
hyperbolic spatial cross-section. Note that the velocity of light c has been taken to be unity. 



As was mentioned earlier, we seek the static, spherically symmetric solutions to the equations ( |2.2| ) and (2.3). To 
this end we assume that the scalar field is the function of r only, i.e. ip — f(r) and the electromagnetic field possesses 
only one component Fio = dAo/dr — A'. 



Under the assumption made above, the solution to the equation (2.3) reads 



F U1 = qP{fp) 



yl — kr 2 



(3.2) 



where q is the constant of integration and P((p) = Putting (3.2) into ( |2.2| ) for the scalar field we obtain the 

equation with "induced nonlincarity" [|L2 J 



(1 - kr 2 ) V " + 2 3fcr V 



2q 2 



R 2 r* v 



This equation can be written in the form 



with the first integral 



0, q 2 =— 
' ^ 16tt 



dz 2 R 2 v 



oz R 



(3.3) 



(3.4) 



(3.5) 



where we substitute z = ^Jl/r 2 — k. Here Co is the constant of integration, which under the regularity condition of Tq 
at the center turns to be trivial, i.e., Co = 0. Finally we write the solution to the scalar field equation in quadrature 



(3.6) 
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In accordance with (3.2) and fl3.5| ) from ( |2.4| ) we find the density of field energy of the system 

Aq 2 P 



T, 



i? 4 r 4 



and total energy of the material field system 

Ef = J T$y^gd 3 x = -8wq J VPdip 



(3.7) 



(3.8) 



Thus, we see that the energy density Tq and total energy Ef of the configurations obtained do not depend on the 
conventional values of the parameter k. As one sees, to write the scalar ((p) and vector (^4) functions as well as the 
energy density (Tq ) and energy of the material fields (Ef) explicitly, one has to give P(p) in explicit form. Here we 
will give a detailed analysis for some concrete forms of P(p)- Let us choose P(p) in the form 



PM = P cos 2 (^) 



with A being the interaction parameter. Inserting ( |3.9| ) into ( |3.4| ) we get the sin-Gordon type equation |1 

Pip Xq 2 P . , . 

1^2 + -r^r sm ( x v) = o 



dz 



AR 2 



with the solution 



p(z) — — arcsintanh[6(2: + z±)], b = , z\ = const 

A R 



(3.9) 



(3.10) 



(3-11) 



Let us analyze the solution ( 3.11 ). It can be shown that lim p = ir/\ for all k. As one sees, the solution possesses 

r— *0 

meaning only in the region where z = y/l/r 2 — k > 0. It means, i n cas e of k = +1 the configuration confines in the 

interval < r < 1. Then the asymptotic behavior of the solution (|3.11 ) can be written as 



0, k = —1, z\ = —1, r — > oo 
0, k = 0, z\ = 0, r — > oo 
0, k = +l, 2i = 0, r>l 



From ( |3.8| ) we find the total energy of the system Ef = —lGirq^/Po/X. For the choice of P(p) in the form 

P = X(a 2 - p 2 ) 2 



with A being the coupling constant and a being some arbitrary constant, from (B.4) we obtain 

4AaV - 4A^ 3 + = 



The equation ( 3.14 ) can be seen as an MKdV one. Indeed, a KdV equation 



can always be converted to 



du „du „d 3 u 

— + au p — +(3t-t -0, 

at ox ox 6 



U p +1 d 2 u n 

Du + a + /3— = 

p + 1 dz 2 



(3.12) 



(3.13) 



(3-14) 



(3.15) 



(3.16) 



if one looks for stationary solution of the form u — u(z) where z — x — Dt. In our particular case p = 2 and the 
equation (3.14) is an MKdV one. The scalar field function in this case has the form 



<p(z) — atanh[VAafe(z + Z2)], 



2q 
R 



(3.17) 
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Taking into account that z = yl/r 
different value of k we get 



k one sees that at the origin lim 02 



a whereas at the asymptotic region for 



0, k — — 1, Z2 — — 1, r — > 00 
0, fc = 0, Z2 — 0, r — > 00 
0, fc = +1, z 2 = 0, r > 1 



(3.18) 



From (3.8) we find the total energy of the system to be = — 167T(7\/Aa 3 /3. A specifi c ty pe of solution to the 
nonlinear field equations in flat space-time was obtained in a series of interesting articles These solutions are 

known as droplet-like solutions or simply droplets. The distinguishing property of these solutions is the availability 
of some sharp boundary defining the space domain in which the material field happens to be located, i.e., the field is 
zero beyond this area. It was found that the solutions mentioned exist in field theory with specific interactions that 
can be considered as an effective one, generated by initial interactions of unknown origin. In contrast to the widely 
known soliton-like solutions, with field functions and energy density asymptotically tending to zero at spatial infinity, 
the solutions in question vanish at a finite distance from the center of the system (in the case of spherical symmetry) 
or from the axis (in the case of cylindrical symmetry). Thus, there exists a sphere or cylinder with critical radius ro 
outside of which the fields disappear. Therefore the field configurations have a droplet-like structure Jl^JT^Jl^ | . 

To obtain the droplet-like configuration we choose a very specific type of interaction function P(ip) which has the 
form @ [cf. FIG. 1] 



where J — Xip, a = 2n + 1, n = 1, 2 • • •. Putting ( 3.19 ) into ( |3.q ) one gets 



|J 2 /--l|=exp[±^(z-z ) 



(3.19) 



(3.20) 



Let us consider the case when | J 2 ! — 1| = 1 — J 2 / ff . Taking the sign in exponent to be "minus" one from ( |3.2C| ) we 
obtain 



tp(z) 



1 1 4qA / v 

l-exp( - — {z-zq) 



a/2 



(3.21) 



Recalling that z = •v/l /r 2 — k from (3.21) we see that at 



> the scalar field tp takes the value tp(0) — ► 1/A and at 
r — > r c = 1/ \J~Zq~- \-k, the scalar field function becomes trivial, i.e., tp(r c ) — ► 0. It is obvious that for r > r c the value 
of the square bracket turns out to be negative and tp(r) becomes imaginary, since a is an odd number. Since we are 
interested in real ip only, without loss of generality we may assume the value of ip to be zero for r > r c , the matching 
at r = r c (i.e., z = zq) being smooth [cf. FIG. 2]. Note that, for k = +1, the scalar field is confined in the region 
< r < 1, as it was in previous two cases. The total energy of the droplet we obtain from (|3.gD has the form 



E f 



X(a-l) 



(3.22) 



As is seen from (3.22), the value of the total energy does not depend on the size of the droplet, it means droplets of 
different linear size share the same total energy. Let us now back to (3.2C) again and consider the sign in exponent 
to be "plus". In this case we find 



<p{z) = - 



I AqX ( \ 



7/2 



(3.23) 



Contrary to the droplet this configuration possesses trivial value up to r = r c = \j \J z 2 + fc, then begins to increase 
taking maximum value at spatial infinity: 

^[l-cxp( — )] <1 



Note that the total energy of the anti-droplet ( [3.23| ) is equal to that of the droplet. Finally we consider another very 
interesting case putting 
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P( (p ) = ^Q 2 -^[A 2 -4Q 2 ^ 



where Q — Rip/2q, a = 2n + 1, n = 1, 2 • • •. Putting (3.24) into (3.6) one gets 



<P(z) = J^[( z - a )( b - z ) 



ict/2 



(3.24) 



(3.25) 



where we put A = b — a and 2zq = b + a. As one sees this configuration is completely localized in the interval 
r G + 1/Va 2 + k) having trivial value out of it [cf. FIG. 3]. The total energy of this configuration is trivial. 



B. Stability problem 

To study the stability of the configurations obtained we write the linearized equations for the radial perturbations 
of scalar field assuming that 



<p(r,t) = (p(r)+£(r,t), £<<^> 



Putting ( 3.26 ) into (2.2) in view of (3^) we get the equation for £(r, t) 

" R £ l-kr\„ 2 - 3kr 2 d t q 2 P^ = Q 



€ + 3 ^" R 2 



R 4 r 4 



(3.26) 



(3.27) 



The second term in (3.27) is zero since we assume the FRW space-time to be static one putting R = constant. 
Assuming that 



£(7-, t) « v(r)exv(-iQ.t), Vt = oj/R 



from ( p.27| ) we obtain 



The substitution 



(1 - kr 2 )v" 



n 2 P 

2 H 1 ipip 



R 2 r 4 



v = 



(3.28) 



(3.29) 



77(C) = r ■ v(r), ( = —=&rcsm(Vkr) 



leads the equation (3.29) to the Liouville one pi 



r }cc + (w 2 -V(tp))v = 0, V(<p) = -k + 



( Vk( y 



R 2 C 4 W(V*C) 



For the interaction term (3^9) we see that 



V{<p)>0, if and only if tanh 2 6(z + z x ) > |^ F^f^)' + \ 



(3.30) 



(3.31) 



(3.32) 



Thus we find that the equations with trigonometric nonlinearity contain stable solutions. Given P(tp) in the form 
( 3.13 ) we find that 



rr/ s n -r ,2w x k( 4 R 2 / sm(Vk() \ 4 a 2 



(3.33) 



As is seen from (3.33) the equations with polynomial type of nonlineari ty to o contain some stable solutions. For the 
droplet-like configurations, i.e., for the interacting term P((p) given by ( 3.19| ), it can be shown that the potential 



lim V{ip) — > +00, lim V{ip) — ► +c 



(3.34) 



beginning with a > 5. It means that the droplet-like configurations ( |3.21 ) with a > 5 are st able f or th e class of 
perturbation, vanishing at r — and r = r c . The same can be concluded for the solutions (3.23) and (3.25). 
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4. CYLINDRICALLY SYMMETRIC CONFIGURATIONS 

A. Solutions in Godel Universe 

In the previous section we studied the possibility of formation of regular localized configuration in homogenous and 
isotropic FRW Universe. Let us now continue our study in the homogenous but anisotropic Universe. In particular 
we consider the model proposed by Godel. The linear element of Godel Universe in cylindrical coordinates reads pl| 

1 /8 
ds 2 = dt 2 -dp 2 + — [sinh 4 Op - sinh 2 fipl<i</> 2 - ^smh 2 Qpd(j)dt - dz 2 (4.1) 

where the constant f2 is related with the angular velocity lu: uj = \/2D,. This form of linear element of the four- 
dimensional homogenous space S directly exhibits its rotational symmetry, since the g^ u do not depend on <fi. It is 
easy to find 

lim J — q = lim — — sinh(2£!o) — > p 

i.e. at uj — > Godel Universe transfers to the flat one. 

As in spherically symmetric case, here too we seek the static solutions to the equations ( |2.2| ) and ( |2.3| ) assuming 
the scalar field to be the function of p only, i.e., tp — ip(p) and the electromagnetic possesses only one component 
Fio = dAo/dp = A' . For the electromagnetic field in this case we find 

F Q1 = 2ftZ)P/sinh(2ftp), D = const (4.2) 



The scalar field equation (2.2) with regards to ( |4.2| ) reads 

g^ + 2ftc th(20p)gg= 8q Y P \ , q 2 =D*/mn (4.3) 
dp 2 ''dp sinh 2 (2fip) ' K ' 



Putting y = ^lntanh(fJp) from (13) one gets 

d 2 <p 



dy 2 1 

with the first integral 



n z P v = (4.4) 



-f = ±4qn^/P + D (4.5) 
dy 

Here Do is the constant of integration, which under the regularity condition of Tq at the center turns to be trivial, 
i.e., Dq = 0. Finally we write the solution to the scalar field equation in quadrature 



J % = ^{y - y ) 



(4.6) 



In accordance with (4.2) and (4.5) from (2T) we find the density of field energy and the total energy of the system 



" sinh 2 (2ftp) 



Ef = 8nq J VPdip (4.8) 

As in the previous case, to write the scalar (tp) and vector (A) functions as well as the energy density (Tq ) and energy 
of the material fields (Ef) explicitly, one has to give P((p) in explicit form. Here ag ain w e will thoroughly study the 
solutions obtained for different concrete forms of P(<p)- Choosing P(<p) in the form (|3.9|), i.e., 



P(^=P cos 2 (^) (4.9) 
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with A being the interaction parameter, from (4.4) we get the sin-Gordon type equation 



o2 

+ 4\q 2 n 2 P sin(\<p) 



The solution to this equation can be written in the form 

/ tanhilpA ' 



<f(p) 



arctan 



(4.10) 



(4.11) 



where po is the constant of integration, giving the size of the system. Without losing the generality we can choose 
a > 0. Then one finds lim p = —ir/\. For p > the field p steadily increases up to 7r/A. In particular, at spatial 

infinity we get lim p> = 0. In this case P(^oo) = 1 which corresponds to the exclusion of interaction at spatial 

p — >oo 

infinity. The total energy of the s ystem in this case coincides with that of in FRW Universe, i.e., Ef — — IQnq^fP^/ \. 
Choosing P(p) in the form ( 3.13| ), i.e., 

(4.12) 



P = \(a 2 - ip 2 ) 2 



from (4.4) we as in FRW case again obtain MKdV type equation with the solution 

. tanh"(f!p) - 1 /- 

pip) — a , nlr ^ , -, a = 4aqV A 

tanh Q (ftp) + 1 



(4.13) 



-a and lim p 

p — >oo 



From ( [4.13| ) it is clear that limp 

as it was in FRW case. The choice of the interaction term P(p) in the form fl3.19| ), i.e., 

%) = J 2 - 4 /»(l-J 2 /f 
where J = A^, a = 2n + 1, n = 1, 2 • • •, leads to the following expression for scalar field 



From (4.8) we find the total energy of the system to be 



(4.14) 



1 - 



/ tanh(ilp) 
Vtanh(r2po) 



-i cr/2 



a = ± 



4A<? 



(4.15) 



where po is an arbitrary constant. For a > the solution possesses physical meaning at p < po and becomes 
meaningless at p > po- Outside the cylinder p = p$ one can put ip = 0. This trivial solution is stitched with the 
solution at p = po under condition p'(po) = 0, which fulfills if and only if i\\q\ > a. Consequently, at p > po the 
Lagrangian becomes physically meaningless, however its limiting value at p — > po — is equal to zero. Continuing it 
at p > po, one can consider the field be totally trivial in this area. Thus we get the droplet-like configuration. The 
field p steadily decreases from p(0) = 1/A to p(po) = with ^'(0) = (for a > 3) and p>'{po) = 0. The total energy 
of the "droplet" is defined as 



E f = 



4irq 



A(a-l) 



(4.16) 



which remains unaltered even in flat space-time (fl = 0). It means that the "droplet" does not feel Godel gravitational 
field. On the other hand, for a < we obtain the configuration (anti-droplet) possessing trivial value at p < po. 
Starting from p = po it begins to increase and at spatial infinity takes maximum value: p(oo) < 1/A. Finally we 
consider the case providing hat-like configuration. Putting 



P{p) 



T 



2/a 



where Q = R/4qfl, a = 2n + 1, 



1, 2 ■ • •. Inserting (4.17) into ( |4.6| ) one gets 
p{y) = Aqn[{y-a){b-y)Y /2 



(4.17) 



(4.18) 



where we put A = b — a and 2yo = b + a. As one sees the configuration is completely localized in the interval y € (a, b) 
with trivial total energy. 
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B. Stability problem 



Let us now study the stability of the configuration obtained. In doing so we consider the perturbed scalar field 
Sip = x{Pi t) that leaves the cylindrical-symmetry of the system unbroken. The linearized equation for the perturbed 
scalar field looks 



1 d ( 



-g dx^ 



2(-.9) 



(4.19) 



Since, for the case considered, <? 00 — (1 — smh 2 flp)/ cos h 2 flp , it is clear that the type of equation ( 4.19[ ) changes on 
the surface where sinhfip = 1. The Cauchi problem for (4.19) is incorrect by Hadamard in the region, where g Q0 < 0. 
In connection with this only the droplet-like and hat-like solutions can be stable by Lyapunov, if they are located in 
the region where g 00 > [Of. Assuming that 



\{p,t) = v(p)exp(-iet) 



from (4.19) we get 



The substitution 



where 



2ftcoth(2ftpV 



1 



sinh flp 2 



2ti 2 q 



1„2 



cosh ftp 



sinli (2f2p) 



P vv ]v = 



vie) = ap) ■ v(p) 



(4.20) 



(4.21) 



(4.22) 



c 



yl — sinh 2 Qp 



cosh Vtp 

leads the equation for perturbed field to the Liouville one 



dp, £ = [4sinh 2 Qp(l - sinh 2 ftp)] 



1/4 



+ [- 2 -^=° 



(4.23) 



Here the effective potential V(ip) takes the form 

n 2 r q 2 P w 



■I- 



(4sinh 6 57/9 — 16sinh 4 fip + 3sinh 2 f2p — 1) cosh 2 f2p ■ 



4sinh 2 fip (1 - sinh 2 Slp) 



(1 - sinh 2 ftp) 1 2sinh 2 ftp 
For the interaction function, chosen in the form ( 3.19 ), one finds 

, 2ct 2 - 12cr + 16 4cr 2 - 12cr + 8 



pip 



Taking into account that lim J 



1 and lim J 

P^PO 



.2J4/ CT a 2j2/* 

— > 0, it can be shown that 



lim V(ip) 
lim VUp) 

P^Po 



-oo for a < 4gA 
-oo for a > 5 



Here we used the fact that g 00 > 0, i.e., sinhSlp < |1|. Thus, as in the previous case we find that the droplet like 
configurations are stable for 5 < a < 4|g|A. Analogically, we can conclude that the hat-like configurations are also 
stable beginning with a > 5. Here it should be emphasized that, contrary to the droplets (anti-droplets) where a may 
be big enough, in case of hat-like solutions a should be small. Otherwise ip itself may be large. Nevertheless, one can 
always choose P in such a way, e.g., 



P = (l/4)((7Q) 2 - 4 / ff [A 2 ~ A(aQ) 2 ^] 



that ip takes reasonable value. 
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5. CONCLUSION 



We obtained the regular particle-like solutions to the scalar field equations with induced nonlinearity in external 
gravitational fields described by Freedman-Robertson- Walker and Godel Universes respectively. Beside the usual 
solitons, a special type of regular localized configurations, known as droplets, have been obtained. It has been shown 
that the droplet-like configurations possess limited energy density and finite total energy and the droplets of different 
linear sizes up to the soliton share one and the same total energy. It is noteworthy to notice that in the spherically 
symmetrical case (i.e., in the FRW Universe) at r c — > oo for k — droplet transfers to usual solitonian solution, 
while for k = ±1 this is not the case. It has also been shown that in FRW space-time equations with different type 
of nonlincarities may contain stable solutions, whereas in case of Godel Universe only the droplet-like and hat-like 
configurations may be stable. It is noteworthy to remark that in FRW space-time with k = +1 the field function is 
confined in the region < r < 1 independent to the choice of interaction function P(<p). Acknowledgements The 
author thanks Prof. G.N. Shikin for useful consultations. 
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Solitons of Induced Scalar Field and their Stability 



FIG. 1. Perspective view of interaction function P((p) providing (anti-)droplet like configurations for different values of a. 
Here (and later on) the thick-line, dash-line, dash-dot-line and thin-line correspond to the value of a = 3, 5, 7, 9 respectively. 
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FIG. 2. Perspective view of droplet-like configurations 
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FIG. 3. Perspective view of the hat-like configurations 




